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Optimizing Traffic Flow 
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We present an economics-based method for deciding the optimal rates at which vehicles are allowed 
to enter a highway. The method exploits the naturally occuring fluctuations of traffic flow and is 
flexible enough to adapt in real time to the transient flow characteristics of road traffic. Simulations 
based on realistic parameter values show that this strategy is feasible for naturally occurring traffic, 
and that even far from optimality, injection policies can improve traffic flow. Our results also allow a 
better understanding of the high flows observed in "synchronized" congested traffic close to ramps. 
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The advent of powerful traffic simulators [|f]-[| has led 
to a spate of new discoveries in the area of vehicular traf- 
fic that agree well with empirical observations From 
moving traffic jams jjj to transitions to states of coherent 
motion § these new results offer insights into a compli- 
cated and socially relevant many-body problem, while 
also suggesting ways of designing controls that can max- 
imize the flow of vehicles in cities and highways ||D • 

The recent discovery of a new state of congested high- 
way traffic, called "synchronized" traffic 0, has gener- 
ated a strong interest in the rich spectrum of phenom- 
ena occuring close to on-ramps . In this connection, 
a particularly relevant problem is that of choosing an 
optimal injection strategy of vehicles into the highway. 
Similar questions have recently been raised regarding the 
most efficient use of the Internet in light of its bursty 
congestion patterns [JlO|. While there exist a number of 
heuristic approaches to optimizing vehicle injection into 
freeways by on-ramp controls, the results are still not 
satisfactory. What is needed is a strategy that is flexible 
enough to adapt in real time to the transient flow char- 
acteristics of road traffic while leading to minimal travel 
times for all vehicles on the highway. 

This paper presents a solution to this problem that 
explicitly exploits the naturally occuring fluctuations of 
traffic flow in order to enter the freeway at optimal times. 
This method, which leads to a more homogeneous traffic 
flow and a reduction of inefficient stop-and-go motions, 
is less sensitive to failures in the control mechanism than 
traditional approaches. 

The basic performance criterion behind this optimiza- 
tion technique is the travel time distribution of vehicles 
jllj which is to be contrasted with the local velocity dis- 
tribution applied in other approaches. The travel time 
distribution is a global measure of the overall dynamics 
on the whole freeway stretch. It allows the evaluation of 
both the expected arrival time of vehicles at a destination 
and its variance, the latter characterizing the likelihood 
that a particular vehicle will have a travel time different 
from the expected one. Both the average and the vari- 
ance of travel times are influenced by the inflow of vehi- 
cles entering the freeway over an on-ramp. From these 



two quantities one can construct a relation between the 
average payoff (the negative mean value of travel times) 
and the risk (their variance) , as is considered in many op- 
timization problems in economics. The optimal strategy 
will then correpond to the point in the curve that yields 
the lowest risk at a high average payoff. In the follow- 
ing, we will show that the variance of travel times has a 
minimum for on-ramp flows that are different from zero, 
implying that traffic flow can be optimized by choosing 
the appropriate vehicle injection rate into the freeway 

In order to obtain the travel time distribution of ve- 
hicles on a highway, we simulated two-lane traffic flow 
via a discretized follow-the leader model, which describes 
the empirical known features of traffic flows quite well 
In our experiments, we extended the simulation to 
several lanes with lane-changing maneuvers and differ- 
ent vehicle types (cars and trucks). We determined the 
travel times of all vehicles by storing the times at which 
they pass two successive cross sections of the road. Em- 
pirical data of this kind could be collected at the gates 
of a toll road. The time difference between the passing 
times corresponds to the travel time needed to traverse 
the freeway stretch between the two cross sections. In 
our simulations, we chose a stretch of length L = 10 km. 

The model distinguishes / neighboring lanes i G 
{!,...,/} of a unidirectional freeway. All lanes are 
subdivided into sites z 6 {1,2, . . . , Z} of equal length 
Ax = 2.5 m. Each site is either empty or occupied, 
the latter case representing the back of a vehicle of 
type a € {1,...,A} with velocity v = uAx/At. Here, 
u G {0, 1, ... , u" lax } is the number of sites that the ve- 
hicle moves per update step At. We have distinguished 
cars (a = 1) and trucks (a = 2). These are characterized 
by different optimal velocities U a (d) with which the ve- 
hicles would like to drive at a distance d to the vehicle 
in front. Their lengths l a correspond to the maximum 
distances satisfying U a (l a ) = 0. The positions z(T), ve- 
locities u(T), and lanes i(T) of all vehicles are updated 
Q every time step At = Is at times T G {1,2,...} 
according to the following successive steps || : 

1. Determine the potential velocities Uj(T + 1) on the 
present and the neighboring lanes i(T) + j with j G 
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{ — 1,0, +1} according to the acceleration law 

Uj(T + l) = [XUaid^T)) + (1 - A)u(T)J . (1) 

Here, the floor function [x\ is defined by the largest in- 
teger m < x, and dj(T) — [Zj(T) — z(T)] denotes the 
distance to the next vehicle ahead (+) on lane i(T) +j at 
position z~j(T). The above equation describes the typi- 
cal follow-the-leader behavior of driver- vehicle units. De- 
layed by the reaction time At, they tend to move with 
the distance-dependent optimal (safe) velocity U a , but 
the adaptation takes a certain time r = A At because of 
the vehicle's inertia. Good results, in the sense of repli- 
cating highway data, are obtained for A = 0.77. 

2. Change lane to the left, i.e. set 

i(T + l)=i(T) + k (2) 

with k = +1, if the vehicle considered can go fastest 
there, i.e. if 

u+i(T + l)>ti (T + l) and >u_i(T + l). (3) 

This is usually the case if the headway on the left lane is 
greatest. Apart from the validity of this incentive crite- 
rion, we demand two extra safety criteria Jl3| | : First, the 
current vehicle position should be ahead of the expected 
position zT^T + 1) of the following vehicle (— ) on the 
left lane, i.e. 

z(T)>z- 1 (T+l). (4) 

Second, the potential velocity on the left lane should not 
be considerably less than the expected velocity uTi (T+l) 
of the following vehicle, i.e. 

ufcCr + ij^gu+tCr + i). (5) 

Once again, realistic results are obtained for q = 0.7. A 
value q < 1 implies that drivers are ready to accept a 
braking maneuver of the follower on the destination lane 
at the next update step. Therefore, the values of q are a 
measure of how relentless drivers are in overtaking. 

Assuming symmetrical ("American") lane changing 
rules for simplicity, a change to the right lane (k = —1) 
is carried out, if the incentive criterion u_i(T + 1) > 
uq(T +1) and > u + i(T + 1) as well as the safety criteria 
z{T) > zZi(T + 1) and u k {T + 1) > quZ^T + 1) are 
fulfilled. Otherwise the vehicle stays on the same lane 
(k = 0). 

3. If the potential velocity Uk(T + 1) on the new lane 
i(T + 1) is positive, diminish it by 1 with probability 
p = 0.001, which accounts for delayed adaptation due 
to reduced attention of the driver and the variation of 
vehicle velocities: 

il with probability p 
ifu fc (T+l)>0 (6) 



The resulting value defines the updated velocity u(T+l). 

4. Update the vehicle position according to the equa- 
tion of motion 

z{T + 1) = z(T) + u(T + 1) . (7) 

Our simulations were carried out for a circular two- 
lane road. After the overall density was selected, vehi- 
cles were homogeneously distributed over the road at the 
beginning, with the same densities on both lanes. The 
experiments started with uniform distances among the 
vehicles and their associated desired velocities. The ve- 
hicle type was determined randomly after specifying the 
percentages of cars (90%) and trucks (10%). 

At a given injection rate, vehicles enter the beginning 
of an on-ramp lane (i.e. a third lane) of 1 kilometer length 
with a uniform time headway and their optimal velocity 
related to the density in the destination lane. Our sim- 
ulation stretch consists of a 1 km long 3-lane part and 
the L = 10 km long two-lane stretch on which the travel 
times are measured. No vehicle is allowed to change to 
the on-ramp from the main road. Injected vehicles try 
to change from the on-ramp to the main road as fast as 
possible, i.e. according to the above lane-changing rules, 
but they do not care about the incentive criterion. The 
end of the on-ramp is treated like a resting vehicle, i.e. 
any vehicle that approaches it has to stop, but it changes 
to the destination lane as soon as it finds a sufficiently 
large gap. If the on-ramp is completely occupied by vehi- 
cles waiting to enter the main road, the injected vehicles 
form a queue and enter the on-ramp as soon as possible. 
Injected vehicles that have completed their 10 km trip 
on the two-lane measurement stretch are automatically 
removed from the freeway (which would correspond to 
uncongested off-ramps adjacent to the lanes). Since our 
evaluations started after a transient period of two hours 
and continued until the 1000th injected vehicle finished 
its trip, the results are largely independent of the initial 
conditions. 

If we plot the average of travel times as a function 
of their standard deviation (Fig. [I]), we obtain curves 
parametrized by the injection rate of vehicles into the 
road. Several points are worth noting: 1. With grow- 
ing injection rate Q rmp — l/2"s (n € {2, 3, ... , 10}), the 
travel time increases monotonically. This is because of 
the increased density caused by injection of vehicles into 
the freeway. 2. The average travel time of injected vehi- 
cles is higher, but their standard deviation lower than for 
the vehicles circling on the main road. This is due to the 
fact that vehicle injection produces a higher density on 
the lane adjacent to the on-ramp, which leads to smaller 
velocities. The difference between the travel time distri- 
butions of injected vehicles and those on the main road 
decreases with the length L of the simulated road, since 
lane-changes tend to equilibrate densities between lanes. 

In addition, the standard deviation of travel times has 
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tend to fill existing gaps and thus homogenize traffic flow. 
This minimum is optimal in the sense that there is no 
other value of the injection rate that can produce travel 
with smaller variance. In particular, gap-filling behavior 
mitigates inefficient stop-and-go traffic at medium den- 
sities. Above a density of 45 vehicles per kilometer and 
lane on the main road (without injection), the minimum 
of the travel times' standard deviation occurs for n = 6. 
The reduction of the average travel time by smaller in- 
jection rates is less than the increase of their standard 
deviation. This result suggests that, in order to gen- 
erate predictable and reliable arrival times, one should 
operate traffic at medium injection rates. Lastly, for the 
case of 40 vehicles per kilometer and lane, the minimum 
of the standard deviation of travel times is located at 
n = 5, while for 35 vehicles per kilometer and lane, it 
is at n = 4. Below 30 vehicles per kilometer and lane, 
vehicle injection does not reduce the standard deviation 
of travel times. This is because at these densities homo- 
geneous traffic is stable anyway, so no stop-and-go traffic 
exists and therefore no large gaps that can be filled ||. 

The curves displayed in Fig. |l| correspond to a given 
density p ma m on the main road without injection of ve- 
hicles. The effective density p e ff on the freeway resulting 
from the injection of vehicles can be approximated by 



Pctl = Pu 



IL 



(8) 



where I = 2, L = 10 km. AT in j is the average number 
of injected vehicles present on the main road and can be 
written as 



' T 



77 • 



(9) 



where N tot = 1000 is the total number of injected vehicles 
during the simulation runs, 7i n j is their average travel 
time, and T to t the time interval needed by all iV to t = 
1000 vehicles to complete their trip. We point out that, 
in addition to these measurements, we used two other 
methods of density measurement which yielded similar 
results. 

We also investigated the dependence of the travel time 
characteristics on the resulting effective densities of vehi- 
cles. As Figure ^| shows, vehicle injection can actually re- 
duce the average travel times of the vehicles on the main 
road, while the travel times of injected cars are about 
the same as those of vehicles on the main road without 
injection. This means that, for given p e ff, one can ac- 
tually increase the average velocity Vmain = i/^main of 
vehicles by injecting vehicles at a high rate without af- 
fecting their travel times. This result follows from the 
increased degree of homogeneity caused by entering ve- 
hicles that fill gaps on the main road, which mitigates 
the less efficient stop-and-go traffic. It might actually 
make sense to diverge a certain proportion of vehicles at 



intermittent times and to reinject them later following 
our methodology. 

Finally, Figure || shows the average of the travel times 
for vehicles in the main road as a function of their stan- 
dard deviation. In contrast to Fig. |l|, the curves were 
computed for the resulting effective densities. This time, 
an increase of the injection rate (which corresponds to 
a smaller number of circling vehicles on the main road 
and a greater proportion of injected vehicles) reduces the 
average travel times. Once again, we observe a minimum 
of the standard deviation of travel times at high vehicle 
densities and medium injection rates. 

In the limit of high injection rates, a traffic jam of max- 
imum density /? max builds up at the end of the on-ramp, 
while downstream of it we find the typical density p out 
related to the universal outflow Q ou t from traffic jams 
HA. We conjecture that the resulting structure con- 
sists of a block of density p max and length L\ containing 
Ni = Pmax-^i vehicles, and a block of density p out of 
length (L — L\) containing {N — N\) = p ut{L — L\) ve- 
hicles at a mean density of p e s = N/L. The expected 
travel time Tmain would be 
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where V on t is the typical velocity emerging downstream 
of a traffic jam, and C is the universal dissolution veloc- 
ity of traffic jams 14 i|. Notice that, for high injection 
rates, the average travel time should grow linearly with 
the mean density p c ft, which is consistent with the results 
displayed in Fig. ^. For decreasing injection rates, travel 
times should increase, since the alternation of congested 
and free flow in the resulting stop-and-go traffic implies 
relevant acceleration times 7a CC and deceleration times 
7dec i n total. 

In conclusion, we have presented a strategy for opti- 
mizing traffic on highways in the sense of higher flows 
and more reliable predictions of individual travel times. 
The applied method is economics-based and resorts to 
the establishment of average payoff versus risk curves. 
Here, the average payoff corresponds to the negative 
mean value of travel times and the risk to their variance. 
The strategy exploits the naturally occuring fluctuations 
of traffic flow in order to allow the entry of new vehicles 
to the freeway at optimal times. Simulations based on re- 
alistic parameter values show that this strategy is feasible 
for naturally occurring traffic, and that even far from op- 
timality, injection policies can improve traffic flow. The 
latter allows an interpretation of the high flows observed 
in "synchronized" congested traffic close to ramps j?J , in 
contrast to the lower average flows in stop-and-go traffic. 

Acknowledgments: D.H. wants to thank the DFG for 
financial support (Heisenberg scholarship He 2789/1-1). 



Hubcrman/Helbing: Optimizing Traffic Flow 



4 



Injection-Rate Dependence of Travel Time for 90% Cars and 10% Trucks 
35 



[1] K. Nagel and M. Schreckenberg J. Phys. I France 2, 2221 
(1992). 

[2] D. Helbing, Phys. Rev. Lett., submitted (1998). 

[3] D. Helbing and B. A. Huberman, Nature, in print (1998). 

[4] M. Schreckenberg and D. E. Wolf (eds.) Traffic and Gran- 
ular Flow '97 (Springer, Singapore, 1998). 

[5] B. Faieta and B. A. Huberman, "Firefly: A synchroniza- 
tion strategy for urban traffic control" (Xerox PARC In- 
ternal Report, 1993). 

[6] P. H. L. Bovy (ed.) Motorway Traffic Flow Analysis. 
(Delft University Press, Delft, 1998). 

[7] B. S. Kerner and H. Rehborn, Phys. Rev. Lett. 79, 4030 

(1997) and Phys. Rev. E 53, R4275 (1996). 

[8] H. Y. Lee, H.-W. Lee, and D. Kim, Phys. Rev. Lett. 81, 
1130 (1998). 

[9] D. Helbing and M. Treiber, Phys. Rev. Lett., in print 

(1998) . 

[10] B. A. Huberman, R. M. Lukose and T. Hogg, Science 

275, 51 (1997); R. M. Lukose and B. A. Huberman, Net- 

nomics, in print (1998). 
[11] K. Nagel and S. Rasmussen, in Artificial Life IV, edited 

by R. A. Brooks and P. Maes (MIT Press, Cambridge, 

MA, 1994). 

[12] We applied sequential update in driving direction, which 
avoids conflicts of vehicles that like to change to the same 
lane and position from both neighboring lanes at the 
same time. For two- lane roads the sequential update leads 
to almost identical results as a parallel update, which is 
usually applied in cellular or lattice gas automata for flow 
simulations and gives realistic results. 

[13] K. Nagel et al, Phys. Rev. E 58, 1425 (1998). 

[14] B. S. Kerner and H. Rehborn, Phys. Rev. E 53, R1297 
(1996). 



Injection-Dependent Travel Times for Different Densities on Main Road 





35 r 




30 ■ 


1 




rimes 


25 ■ 


•avel 


20 ■ 


t— 




50 
n: 


15 - 






> 

< 


10 ■ 

5 L 



Main Road — 




Injected Cars 






,-''50 veh/(km lane^,,-- 




45 veh/(km lane) 


.-*,- yteC,.. 






40 veh/(km lane) 


*s-fe-i:. ; " 35 


veh/(km lane) 


*^*£.',-—30 veh/(km lane) 


<sV- 25 veh/(km lane) 


20 veh/(km lane] 





012345678 
Standard Deviation of Travel Times (min) 

FIG. 1. Average and standard deviation of the travel times 
of vehicles on the main road and from the ramp as a function 
of the injection rate for various vehicle densities on the main 
road (measured without injection). With increasing injection 
rate, the average travel times are increasing due to the higher 
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FIG. 2. Average travel times of vehicles on the main road 
as a function of the resulting effective vehicle density on the 
freeway, for various injection rates. In the limit of high in- 
jection rates, one observes the predicted linear dependence 
of average travel times on effective density, see Eq. ([ji^). In 
contrast to the vehicles on the main road, the travel times of 
injected vehicles did not depend on the injection rate. How- 
ever, when we checked what happens if the vehicles on the 
main road try to change to the left lane along the on-ramp 
in order to give way to entering vehicles, we found that both, 
injected vehicles and the vehicles on the main road, profited. 
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FIG. 3. As figure [lj but as a function of the resulting ef- 
fective density on the freeway. We find shorter travel times 
at high injection rates because of the homogenization of traf- 
fic. The standard deviation of travel times is varying stronger 
than the average travel time, which indicates that medium in- 
jection rates are the optimal choice at high vehicle densities. 



